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OBJECTIVES 

The objective of our research effort is to develop a knowledge base in computational electromagnetics that 
can impact a wide variety of technologies. To achieve this objective we will perform basic research on 
two fronts: 

• Fast integral equation solvers for analyzing scattering and propagation problems in a numerically exact 
manner with reduced computational complexity and memory requirements. 

• Accurate partial differential equation solvers that minimize grid-dispersion errors, and that are devoid 
of spurious solutions. 

In addition, we will augment the above by hybridizing differential equation and integral equation solvers to 
exploit their respective strengths. This will allow us to solve scattering problems involving complex 
environments like a scatterer buried in an inhomogeneous, anisotropic media. 

In order to seek the solution of extremely large or quasi-optical problems, we will hybridize the 
aforementioned solvers with high-frequency methods such as the shooting and bouncing ray (SBR) method 
that is used in XPATCH. 

The multiple scattering T matrix approach has been used to analyze scattering from a large number of 
randomly distributed particles. We will improve on this algorithm for particles of complex shapes and use 



the resulting codes to calculate parameters for the elements of Mueller matrices to be used in transport 
equation like the radiative transfer equation for simulation of wave propagation through random media. 

Finally, to harness maximally the throughput of modem computers, we will seek the parallelization of the 
aforementioned algorithms on massively parallel processor (MPP) machines. This will allow us to solve 
problems of unprecedented sizes and complexities. 

At UT Austin, the scientific objective of this research unit is to develop postprocessing algorithms that are 
tightly coupled to CEM engines to enhance phenomenology understanding and provide linkage to the end 
user community. In particular, we are addressing postprocessing issues in both radar signature prediction 
and antenna design for complex platforms. In the signature prediction area, we are developing 
extrapolation and interpolation algorithms in frequency and angle to accelerate the generation of range 
profiles and inverse synthetic aperture radar (ISAR) imagery using the Fast Illinois Solver Code (FISC). In 
the antenna design area, we are developing the antenna synthetic aperture radar (ASAR) concept to 
pinpoint locations of dominant antenna-platform interaction from CEM data. 


STATUS OF EFFORT/ACCOMPLISHMENTS/NEW FINDINGS 

On the frequency domain fast solvers front, we have: 

• We have continued our focus on improving the Fast Illinois Solve Code (FISC) which has a user base 
of over 400 distributed by SAIC/DEMACO. To stress the code for large scale computing, we put 
FISC through a number of tour de force computations. Recently, we have put FISC through a broad 
bandwidth calculation where we can simulate the time-domain scattering from a full-size aircraft with 
frequency content up to 1 GHz. 

• We continued to make progress with the parallelization of MLFMA with the development of ScaleME. 
By the end of April of2001, ScaleME was able to solve a 10 million unknown problem. The code is 
about 7-9 times faster that FISC (Fast Illinois Solver Code) that solved a similar size problem about 
over two year ago. The speedup compared to FISC essentially comes from the good scaling property 
of ScaleME for massively parallel computing. We used 128 processors on the NCSA SGI 02K. 

• In addition, we have looked at ways to solve Maxwell’s equations with fast solvers all the way from 
static to microwave frequencies. This simulation capability will encompass circuit theory as a special 
case, allowing antenna impedances and circuit phenomena to be more accurately modeled. 

• We have also developed a new fast solver when the scatterer is placed on top, or embedded in a 
layered earth. This capability will have significance in studying the radar returns from targets placed 
under tree canopy, or in studying scattering from buried objects. With this code, we can solve 
problems with over a million unknowns easily. 

In the time domain and genetic algorithm fronts, the progress are: 

• A new PWTD kernel suitable for the analysis of low frequency problems was developed. Just like all 
all fast multipole like schemes, the classical PWTD kernel suffers from a "low frequency breakdown 
problem". This new kernel effectively resolves this issue wiityhin the context of PWTD schemes. 

• A new low frequency PWTD enhanced time domain integral equation solver was developed. This new 
solver relies on loop tree decompositions and permits the analysis of transients on electrically small 
structures using very large CFL numbers. 

• A new PWTD enhanced time domain integral equation solver for lossy penetrable media was 
developed. When hybridized with our previous PWTD solvers for perfectly conducting scatterers this 
new solver will permit the analysis of scattering from coated structures. 

• Inroads were made in the development of PWTD schemes for analyzing layered media problems. 

• New FFT-based time domain integral equation solvers were developed. These solvers constitute 
extensions of the adaptive integral equation scheme to the time domain. While these solvers are not 
expected to compete with PWTD solvers when applied to the analysis f arbitrarily shaped structures, 



they do compete favorably with PWTD technology when the scatterer is a penetrable volume and/or 
quasi-planar. 

On higher-order methods for computational electromagnetics, we have: 

• Developed a time-domain finite element method (TDFEM) for scattering and radiation analysis. In this 
work, we have (1) developed a set of novel orthogonal vector basis functions which yield a diagonal 
mass matrix thus avoiding a matrix inversion in each time step, (2) developed a general approach to 
investigate the stability of the TDFEM approaches and estimate the maximum time step for 
conditionally stable TDFEM schemes, (3) developed an approach to deal with dispersive media which 
enables wideband simulations, (4) developed a perfectly matched layers (PML) for the TDFEM 
simulation of open-region problems, and (5) developed a hybrid approach that combines the TDFEM 
and time-domain integral equation (TDIE) methods for scattering and radiation analysis. 

• Developed a novel hybridization for the high-order finite element method and a high-order multilevel 
fast multipole algorithm (MLFMA) for scattering by complex targets with complex material 
compositions. The hybrid method starts from an absorbing boundary condition and iteratively 
constructs an adaptive absorbing boundary condition that is numerically exact for a given problem. 
Since it solves only a purely sparse finite element matrix, its efficiency is at least an order of 
magnitude higher than traditional finite element—boundary integral methods. 

• Developed a novel physics-based preconditioner for the finite element-boundary integral method. 
Instead of using a mathematical preconditioner, this method construct a preconditioner using an 
appropriate absorbing boundary condition. By using a Krylov subspace-based iterative slover, it can 
achieve a speedup of more than two orders of magnitude over the traditional finite element-boundary 
integral methods, especially when high-order discretization is used. 

• Developed a high-order multilevel fast multipole algorithm (MLFMA) for scattering by objects that 
consist of both perfect electric conductors and piecewise homogeneous materials. This method is based 
on the PMCHWT formulation for the fields in penetrable regions and the combined field integral 
equation for conducting regions. It permits the simulation of complex targets such as coated objects. 

On the post-processing front done at UT Austin, we have: 

• A preconditioner based on the wavelet basis has been developed to accelerate the convergence of 
iterative solvers for large-scale three-dimensional problems. Numerical results on a 3-D cavity showed 
that the wavelet preconditioner outperform the traditional space-domain block-banded preconditioner. 
The computational cost of the preconditioner can be kept under 0(NlogN ), making it compatible with 
the multi-level fast multipole method. 

• The antenna-platform interaction problem has been extended to account for both array mutual coupling 
and platform effects. A sparse model to describe the active element patterns of an array has been 
developed, and an algorithm to extract the model parameters has been devised. Simulation results 
using die latest FISC-radiation code showed that the model can accurately represent very complex 
array responses. Measurements were conducted using a 7-element smart antenna testbed at the 
University of Texas to verify the simulation. 
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Interacted with Raju Numburu at ARL, Adelphi, Maryland. (W. Chew) 
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Fast Solvers for Electromagnetics (W. Chew) 

DEMACO, Champaign, IL 
Dr. D. Andersh 

Worked with Lockheed Martin Corporation to transfer the inlet-simulation technology to industry (J. Jin). 

• Transitions 
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multipole algorithm that can solve complex structures. We call this code FISC (Fast Illinois Solver Code). 
We hope to add new features to the code to enhance its capability. Many copies of the code have been 
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NEW DISCOVERIES, INVENTIONS, OR PATENT DISCLOSURES 

• Developed FISC to a higher level of sophistication with reduced setup time, and capability to solve 
poblems with close to 10 million unknowns as well as performing broadband time domain 
computation. 

• Implemented a highly portable and parallel version of multilevel fast multipole algorithm (MLFMA) 
and tested on SGI Array and Linux cluster solving problem up to 10 million unknowns. 

• Developed a new fast algorithm for layered media using fast inhomogeneous plane wave algorithm. 

• Developed a high-order time-domain finite element method. 

• Developed a novel preconditioner for the finite element-boundary integral method. 

• Developed multilevel plane wave time-domain algorithms for analyzing large scale acoustic scattering 
phenomena, and electromagnetic phenomena involving perfectly conducting as well as penetrable 
targets. 

• Developed fast global boundary conditions for FDTD simulations based on PWTD technology. 

• Developed PWTD based techniques for analyzing nonlinear phenomena. 

• Developed diploidy GAs for controlling adaptive antennas. 

• Developed multiparameter model-order reduction techniques. 

• Developed low-frequency methods both in the time and the frequency domain. 

• Developed higher-order accurate MOM and FEM methods. 

• Developed efficient method for the engine inlet problem. 

• H. Choo and H. Ling, "Microstrip antennas and methods of designing same,” provisional patent 
application, docket no. 5119-09201, submitted December 2001. 
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